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We set forth a polarization-sensitive quantum-optical coherence tomography (PS-QOCT) tech- 
nique that provides axial optical sectioning with polarization-sensitive capabilities. The technique 
provides a means for determining information about the optical path length between isotropic re- 
flecting surfaces, the relative magnitude of the reflectance from each interface, the birefringence of 
the interstitial material, and the orientation of the optical axis of the sample. PS-QOCT is immune 
to sample dispersion and therefore permits measurements to be made at depths greater than those 
accessible via ordinary optical coherence tomography. We also provide a general Jones matrix the- 
ory for analyzing PS-QOCT systems and outline an experimental procedure for carrying out such 
measurements. 

PACS numbers: 42.50.Dv, 42.65.Lm 



I. INTRODUCTION 



Optical coherence tomography (OCT) has become a well-established imaging technique 



in ophthalmology |a| , intravascular measurements 
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with applications 

3, l3l , and dermatology 8] . It is a form of range- finding that makes 
use of the second-order coherence properties of a classical optical source to effectively section a reflective sample with 
a resolution governed by the coherence length of the source. OCT therefore makes use of sources of short coherence 
length (and consequently broad spectrum), such as superluminescent diodes (SLDs) and ultrashort-pulsed lasers. As 
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broad bandwidth sources are developed to improve the resolution of OCT techniques, material dispersion has become 



more pronounced. The deleterious effects of dispersion broadening limits the achievable resolution as has been recently 
emphasized ^ . 

To further improve the sensitivity of OCT, techniques for handling dispersion must be implemented. In the partic- 
ular case of ophthalmologic imaging, one of the most important applications of OCT, the retinal structure is located 

10| . Dispersion increases the width of the coherence 



behind a comparatively large body of dispersive ocular media 

envelope of the probe beam and results in a reduction in axial resolution and fringe visibility . Current techniques 
for depth-dependent dispersion compensation include the use of dispersion-compensating elements in the optical set- 
up 2 2 or employ a posteriori numerical methods ^ 14 1. For these techniques to work, however, the object 
dispersion must be known and well characterized so that the appropriate optical element or numerical algorithm can 
be implemented. 

nil 

Over the past several decades, a number of non-classical (quantum) sources of light have been developed llallal and 
it is natural to inquire whether making use of any of these sources might be advantageous for tornographic i maging . 
An example of such a nonclassical source is spontaneous parametric down-conversion (SPDC) [rUlsHlS. 2ol 2li 2^. 
a nonlinear process that produces entangled beams of light. This source, which is broadband, has been utilized 
to demonstrate a number of interference effects that cannot be observed using traditional classical sources of light. 
We make use of this unique feature in quantum-optical coherence tomography (QOCT) 2^, where fourth-order 
interference is used to provide range measurements analogous to those currently obtained using classical OCT, but 
with the added advantage of even-order dispersion cancellation ^ 2^ 2^. We have recently demonstrated the 
dispersion immunity of these tomographic measurements in comparison to standard optical coherence tomography 
techniques |27|. 

In this paper, we present a method for polarization-sensitive QOCT (PS-QOCT) measurements, where one can 
detect a change in the polarization state of light reflected from a layered sample [28| . This state change arises from 
scattering and birefringence in the sample and is enhanced in specimens that have an organized linear structure. 
Tissue that contains a high content of collagen or other elastin fibers, such as tendons, muscle, nerve, or cartilage, 
are particularly suited to polarization-sensitive measurements ^ . A variation in birefringence can be indicative of a 
change in functionality, integrity, or viability of biological tissue. 
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II. GENERAL MATRIX THEORY FOR PS-QOCT 



We present the theory for PS-QOCT according to the simphfied diagram for an experimental setup given in Fig. ^ 



Using a Jones-matrix formahsm similar to that in Ref. 



3] 



, we start by defining a twin-photon Jones vector Jj, 



(1) 



ai{uj') Gi 

where as(w) and ai(aj') are the annihilation operators for the signal-frequency mode lu and the idler-frequency mode to' , 
respectively. The vectorial polarization information for the signal and idler field modes are contained in e^ , (j = s, i). 
For example, if we utilize coUinear type-II SPDC from a second-order nonlinear crystal (NLC) to generate a pair of 
orthogonally polarized photons, the reduce to the familiar Jones vectors: 



e, = 



(vertical) 
(horizontal) 



(2) 



for signal and idler, respectively. 

The twin photons coUinearly impinge on the input port of a polarizing beam splitter (PBS) from which the signal 
photon is reflected into the sample arm and the idler photon is transmitted into the reference arm. We assume that 
the polarization in each arm is independent of that in the other arm until the final beam splitter. We also assume 
that all optical elements within the interferometer are linear and deterministic. 

The delay accumulated by the signal and idler beams in each path is represented by the 2x2 matrix 



D 



(3) 



di(cj) 
d2{uj') 

where di(ti;) and 6.2(1-0') represent the Jones matrices that describe the delay for the sample and reference arms, 
respectively. The polarization state in each arm is represented by the matrix 



U = 



(4) 



Ui(w) 

U2(CJ') 

In the experimental realization of PS-QOCT, the matrix Ui{uj) represents the properties of the sample plus any 
other polarization elements in the sample arm, whereas U2(aj') represents the user-selected polarization state in the 
reference arm. 
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The mixing of the polarizations from each path, which occurs at the final beam splitter, is represented by the 
transformation matrix 



-BS — 



(5) 



T31 T32 

T41 T42 

where each element T^; (fc = 3, 4 and Z = 1, 2) is a 2 x 2 matrix that represents the mixing of independent polarization 
modes from the input paths ® and @ prior to detection in paths @ and ®. For example, T31 is the Jones matrix that 
represents the transformation of input spatial mode 1 into output spatial mode 3. 

The Jones vector Jout that describes the field operators at the output of the final beam splitter, can be computed 
from the product of the previously defined matrices in Eqs. - ® as 



TbsUDJ, 



(6) 



T3iUidi T32U2d2 
T4iUidi T42U2d2 

^From this equation, the fields in paths @ and ® arriving at each of the two detectors can be written in the time 
domain as 



E^+)(t3)= / dwe-''^*^as(w)e33+ / dw' e^'-^'*^ ai(w') 63: 



(7) 



E^+)(i4)= / dwe-''^*''as(w)e4s+ / dw' e"'"''** ai(w') e4i, (8) 



where 

635 = T31 Ui di Gs e3i = T32 U2 d2 Gi 

e4s = T41 Ui di Bs e4i = T42 U2 d2 e; (9) 

describes each of the transformations of the signal and idler polarizations that contribute to the final fields in @ and 
® at the detectors. 

/^From the field at each of the detectors, the two- photon amplitude can be written as 

AMt3,U) = {0\Ei+\ts)E[t\u)\^), (10) 

where j and k represent two orthogonal polarization bases such as horizontal/ vertical (H/V), right /left circular (R/L), 
or (45°/-45°). The ket j^I') represents the two-photon state at the output of the nonlinear crystal, defined as 

l^-) = [ dn^in)aliujo + n)al{ujo-^)\0), (11) 



where <I>(fi) is the state function 



3^ 
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that governs the spatio-temporal properties of the signal and idler photons at 
angular frequency uq^^- The state function is given by ^{il) = Lsinc[AA:^(il)L/2], where L is the crystal length and 
Akz{Q) is the wave- vector mismatch in the z- or phase-matching direction. If the state function is symmetric 

about the center frequency or $(— fi) — $(ri), the SPDC spectrum becomes |$(f2)p. 

We assume that the detection apparatus is slow and independent of polarization so that the final coincidence rate R 
is computed as the magnitude-square of the two-photon amplitude summed over each polarization mode, integrated 
over time: 

I dh I di4^^|A,fe(t3,i4)|'. (12) 

j k 

III. SIMPLIFIED CONFIGURATION FOR PS-QOCT 

It is now useful to consider a specific experimental configuration from which we can define expressions for the Jones 
matrices in Eqs. © " © and calculate the quantum interferogram. One particular experimental configuration for 
PS-QOCT is based on the Mach-Zchndcr interferometer and is shown in Fig. |2| 

The elements of the Jones matrix in Eq. ijJJ), representing the delay in each path, are given simply by di (w) = I e"^^=/^ 
and d2(a;') = le"^ ^i/c^ where I is the identity matrix, c is the speed of light in the medium, and Zg, are the path 
lengths in the sample and reference arms, respectively. We define a path-delay difference r = (zr — Zs)/c between 
the reference and sample arms that becomes our experimental parameter in the final expression for the measured 
coincidence rate R(t). We assume that the final beam splitter faithfully transmits and refiects each input polarization 
mode. The elements of Tbs in Eq- © thus become T31 = T42 ~ tl and T32 = — T4j^ = r* I, where f designates a 
matrix transpose and conjugation, and t and r represent the amplitude transmittance and reflectance of the beam 
splitter, respectively. 

Without having to specify the elements of U in Eq. Q), an expression for the measured coincidence rate as a 
function of the path-delay difference r is calculated to be 

i?(T) cx Ao - Vbs Re [A(2t)] , (13) 

where Aq and A(t) are defined as 

Ao = I dnmn)\^[ellj\{ujo + n)Vi{ujo + n)e,]{elvlV2e,) (14) 



and 



A(r) = / dn Fiujo + n) F*{ujo - fl) e 



-iOr 



(15) 



representing constant and varying contributions to the quantum interferogram, respectively. The function F{ll}), 
which includes all of the sample properties, is given by 



(16) 



The parameter Vbs = 2(|r'p|tp)/(|r|'* + 1^1"*) in Eq. H13|l represents a visibility factor for a lossless beam splitter with 
arbitrary transmittance (Vbs = 1 when |rp = |ip ~ 1/2). We assume that the optical elements in the reference arm 
are frequency independent across the bandwidth of the light-source spectrum. Equation H15|l is a generalization of 
Eq. (8) in Ref. 23], where the function F contains polarization-dependent information about the sample. 

It is clear from Eq. (|15|l that the sample is simultaneously probed at two frequencies, luq + and ujq — Q, and 
that for a frequency-entangled two-photon state such that produced by SPDC, even-order dispersion from the sample 
is cancelled in PS-QOCT. The effectiveness of even-order dispersion cancellation is related to the spectrum of the 
source used for SPDC. Since we assume a cw-pump source in Eq. (|ll|l . this leads to signal and idler photons that 
are exactly frequency anti-correlated. In this case, even-order dispersion cancellation is perfect. As the bandwidth of 
the SPDC pump source is increased, the requirements for exact frequency anti-correlation are relaxed and dispersion 
cancellation is degraded. It is apparent that the delay r can be adjusted to target specific regions in the sample from 
which polarization information can be extracted by scanning the parameters of the user-selected polarization rotator 
U2. This experimental method is similar to those used in quantum ellipsometry, as discussed in Ref. 

In the following section, we consider a specific construct for Eq. Q that defines the optics in the experimental 
setup represented in Fig. |5] Once we derive an expression valid for an arbitrary sample, we consider several special 
cases in an effort to understand the nature of the information contained in the quantum interferogram. 



[iza t 
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IV. ROLE OF POLARIZATION IN PS-QOCT 



To facilitate the description of PS-QOCT, we make use of the Pauli spin matrices 



CTi = 



1 

1 



0-2 



-i 

1 



0-3 = 



1 
-1 



(17) 



from which any 2x2 Hermitian matrix can be defined as A — cqI + c ■ cr, where c = (ci, C2, C3), cr — (cri, cr2, (T3), 
and c • cr denotes the scalar product of vectors c and cr. We first define Eq. Q for N reflective layers, where each 
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reflection is assumed to be isotropic, then consider the special cases of a single and double reflector. 



A. N reflective layers 

We begin with a sample comprised of N reflective layers, each with an interface defined by a reflectance matrix 
rm(Li;), as shown in Fig.|21 The material properties of each layer are represented by a Jones matrix S„i that is assumed 
to be deterministic. The Jones matrix is a product of: an average phase delay; rotation matrices Rm to account 
for the orientation of the fast axis of the sample with respect to the horizontal axis; and the Jones matrix hm for 



S^- If we i; 



a linear retarder with its fast axis oriented along the horizontal axis |32j . If we ignore losses due to absorption, then 
for a single layer of thickness dm = {zm — z^-i), the Jones matrix is given by 

= e'^"(^'"^")B™(d™,a„,c^) (18) 

where Amidm,^^) = (-uhdm/c is the average phase delay of the signal photon at angular frequency lu attained by 
propagating through a layer with average refractive index n = (tIq + «c)/2. The single-pass retardation in the layer 
is given by Sm{dm, uj) — uj Andm/c where Art = Uo — ric is the difference in refractive indices along the fast and slow 
axes of the medium. The rotation matrix R,„(am) and the Jones matrix b„j((5,„) for the linear retarder are given by 

R™(a™)=e-'"-'^^ and b„((5„) = e'f^-/^) <T3 ^ (19) 

respectively, where in general e"''*'"' = (cos 7) I — (isin7) cr. 

A complete transfer function describing the entire sample in Fig.|31is therefore constructed as 

AT 

H(<j-') = S1S2 • • • S„i_iS,„ r,„ S„jS„i_i • • • S2S1 

m=0 
N 

= ^ B(")r„,B("), (20) 

where (/j*-'^-' — 0, -B^"-* = I, and the tilde operation denotes a matrix that takes an argument at a negative angle, viz. 
Sm(a) ~ S„i(— a). The accumulation of all phases up to interface m is given by 

^^Az(d,,c.) (21) 
1=1 

and the accumulated effect of birefringence up to interface m is expressed via 



B(™) ^[]Bi(d,,a,,w). (22) 



1=1 
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Since the principal axes of the layers are generally unknown, a quarter-wave plate (Q) set at 45° is used to convert 
the signal photon to left circularly polarized light. A quarter- wave plate is a polarization element, so that we must 
also include its Jones matrix in the final expression for Ui(cj) (see Fig.0, 

Ui(w) = Q(45)H(u;)Qt(45) 

N 

= ^ e'2^""'Q(45)B(™)r„,B(")Qt(45), (23) 

m=0 

where the quarter-wave plate at 45° is defined as Q(45) — e'^'''^ and Q^(45) = Q(— 45) — e^'^^'^. 

The reference arm only contains a half- wave plate that can be set at an angle 6 to the horizontal axis (see Fig. ^ , 
so that U2 can be written as 

U2(a;') = R(0) e'^'^^ R^{e), (24) 



where, for example, given a vertical input polarization, 29 = 0° selects vertical polarization and 26 — 90° selects 
horizontal polarization. 

We can now write the expression in Eq. assuming frequency- independent isotropic reflection, i.e. = ""s, 

as 

AT 

F{iv) = ejuj ^ e'2^""' Q(45)B(™)r„B(™)Qt(45)e, 

m=0 

N 

= elvl e'2^""V^u^(a;) (25) 

ni—O 

where Um(w) = Q(45) B^™) 0-3 B^™) Q'l'(45) eg. For the sample provided in Fig. |31 the general Eqs. idJ and ifT^ 
become 

N N „ 

Ao = ^ 5] < dn ei2b<"'(^o+n)-^(")(..o+a)] (^^ + f^) + (26) 

m=On=Q 

and 

N N „ 

A(r) = ^ / rf^^l*(f^)l'e'2['^""'(-«+")-'^'"'(-''-f^)l^:*(cc>o-r!)F„,(^o + f^)e-'"", (27) 

m=On=0 

where Fm,{uj) = ej Um{i^)- By substituting Eqs. (|26|l and (|27(l into Eq. H13(l . we construct the final expression for 
the quantum interferogram. In the following sections, we investigate several special samples to explain the features 
contained in the quantum interferogram and to determine a method for extracting sample information. 
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B. Single reflective layer 

If we consider the special case of a single isotropic reflector buried under a birefringent layer of thickness di = zi, 
Eq. H20|l can be written as 

niuj) = SiriSi 

= (28) 

whereupon Eq. H16|l becomes 

F{uj) — irie'^^^ (g| [(i sin(5 sin2ai) I + (cos(5) cri + (sin(5 cos2q;i) 0-3] e^) 

= ine'^^^FiH, (29) 

with 

Fi (uj) = cos S{lu) cos 261 + smS{uj) sin 29 6^'"^ . (30) 

We have made use of the properties of the Pauli spin matrices and the fact that eg = cr^es and ej = crieg. 
For a single reflector, Eqs. (|26|) and H27|) therefore become 



Ao = \n\' J dnmn)\' (31) 

and 

A(r) = |ri|2 / ^^^7|$(^7)|2e'2['^l("o+")-'^l('^«-f^)li^l(^Jo + ^^)Fl*(^^o-^^)e~'"^ (32) 



respectively. The varying term can be further simplified if we expand the propagation constant P{w) in the expression 
for the phase delay, Ai(ci;) — ujhzi/c = f3{uj) z\. The quantity /3(wo + f^) is expanded to second order in so that 
/3(a;o + r2) « /3o + + ^/3"r2^, where /3' is the average inverse of the group velocities Wq and at wq, and /3" 
represents the average group- velocity dispersion (GVD). It is clear that second-order dispersion is cancelled in the 
simplified expression for the varying term, which is given by 

A(r) = |ri|2 \ d^\^{Sl)'^ F^[ujq^^)FI{ujq~^)q-''^^-'^^^'^^\ (33) 



Figure ^ displays the expected curves for a single reflector buried beneath 120 /im of quartz with = 1.54661, 
o — 1.53773, and |rip — 1, using the scheme shown in Fig. [3 For this simulation, we ignore the frequency dependence 



oi S [6 = 5{u!q)], assume that / dil |$(r2)p = 1, and select the fast axis of the quartz sample to be aligned with the 
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horizontal axis in the laboratory frame so that ai = 0. The SPDC spectrum is calculated explicitly via solutions to the 
phase- matching conditions using published Sellmeier equations for BBO. We are interested in the particular case of 
degenerate, collinear type-II phase matching. The top two curves represent the expected coincidence rate, normalized 
by Aq, when the sample photon is mixed with a vertically polarized (i?v, dash-dot curve) or a horizontally polarized 
(i?.H, dashed curve) photon from the reference arm. The solid curve represents the re- normalized total coincidence 
rate [Rt = (i?v + ^^h ^ Ao)/Ao] from which the reflectance of the layer can be recovered. 

The material properties are revealed by the relative values at the center of the dip where r = — 4/3'zi. At this 
value of T, the path-length difference between the arms of the interferometer is zero and there is maximal quantum 
interference. If we neglect any frequency dependence in the birefringence, we can substitute Eq. (|30|l into Eq. H33() 
and write an expression for the coincidence rate at the center of the dip as 

A(r = -4/3'zi) = = \ cos5 cos 261 -I- sin J sin2e e^'°''\^ . (34) 

In the particular case when we select the linear-rotator angle 29 to be either 0° or 90°, corresponding to a polarization 
of the reference photon that is horizontal or vertical, respectively, we obtain 

Ah = cos^J 

Av = sin^S. (35) 
It is possible to determine the value of 5, or the birefringence An, by forming a ratio of these rates at Az = 0: 

= uq An zi/c. (36) 

We can neglect the frequency dependence of S{u!) = ujAn{u!) zi/c — Sf3{u!)zi when 5'{uJo) Afi <C 1, where Aft is the 
bandwidth of the SPDC spectrum. In this limit, the width of the interference dip is larger than the delay between 
the signal and idler fields resulting from the birefringence of the layer. If the bandwidth of the SPDC spectrum is 
increased, the opposite limit can be realized, namely 6'{uJo) Afl ^ 1. In this case, the interference pattern comprises 
of three regions: the expected central dip at r = — 4/3'zi provides the value of S{uJo) = S a,s in Eq. 136(1 : and two 
additional satellite interference patterns centered at r = —[4/3' ± 2(5/3'] zi, where Sf3' is the coefficient of the first-order 
expansion of Sf3 in fi, and provide information about the group- velocity dispersion in the layer. 

Since we choose the linear-rotator angles 29 to be either 0° or 90° , any dependency of the coincidence rate according 
to the orientation angle ai is lost. It is possible, however, to extract the value of ai by using a technique that is 
analogous to null ellipsometry. In the reference arm, if we combine the linear-rotator used to rotate the linear input 



S = tan ^ 



Av 
Ah 



(37) 
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polarization state ej with a quarter-wave plate to transform the linear polarization into a general elliptical state, it is 
possible to exactly match any polarization state in the sample arm. This transformation is given by 

cos 29 + i cos 2(0 — i 
sin2e' + isin2(0- 6 

where cj) is the angle of the quarter-wave plate and 9 is the angle of the linear rotator fast axes with respect to the 
horizontal axis. In the special case when cj) = 29, we revert to the case of a single linear rotator as in our previous 
example. 

When the polarization in the reference arm is selected by this cascade of polarization elements, we can write Eq. (|30f) 

as 

Fi (w) = {cos 5{uj) [cos 261 - i cos 2((/) - 9)] + sin (5(w) [sin 261 - i sin 2((?!) - 9)] e^'"^ ] . (38) 

If the values of (f) and 9 are adjusted so that the polarization state in the reference arm is exactly orthogonal to that in 
the sample arm, [-Fil^ = and the coincidence count rate will be maximized. The value for a can then be determined 
by solving the following conditions of orthonormality, namely, the real and/or imaginary parts of Eq. H38fl must equal 
zero 

cos(5 cos26' + sin(5 sin26' cos2q;i — sin(5 sin2((/) — 6') sin2Q;i — (39) 
— cos(5 cos2(0 — 6*) -t- sin(5 sin26' sin2ai — sin5 sin2((/) — 0) cos2ai = 0. 

If the value of 5 is known, then only one of these equations is required. 

C. Two reflective layers 

A sample with reflections from two surfaces separated by a birefringent material can be expressed as 

= ro + SinSi 

= ro + e'^'^iBinBi (40) 

where the subscripts and 1 denote the first and second boundaries, respectively. In this case, the function in Eq. I|16|) 
becomes 

F{uj) =iroFo + irie^^^ Fi, (41) 
where Fq = cos 26* and Fi has been provided in Eq. (|30|l . 
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For two reflectors separated by a birefringent medium, the constant and varying contributions from Eqs. (|26(l and 
(|?7jl become 

Ao = |ropy" dn\^{n)f + \rif J d^! + r* n e'^/^o-i J dn\'S?{n)\^ului{LUo + fl) &^'^^'"+>^"^'^'' + cc. (42) 

and 

A(r) = |ro|2 gf") (r) + |ri 5(1) (r - 4/3'zi) + n g^' " 2/3' ^1) e'^''"^^ + cc, (43) 

respectively, where the subscript d denotes a contribution that is subject to even-order dispersion and cc indicates 
the complex conjugate, with 

5^"^^) = J dnmn)fFraiu;o + n)F:^{LUo~n)e-'''^ 

and 

The first two terms in Eq. (|42|l arc contributions to the constant coincidence rate arising from each of the two interfaces 
in the material. The third term introduces a contribution only when these interfaces have a separation that is less than 
the coherence length of the signal photon. The first two terms in Eq. H43|l represent dips arising from reflections from 
the first and second surfaces. The third term, which appears midway between these dips, arises from the interference 
between probability amplitudes associated with each of these refiections. 

Figure |S1 provides numerical results for a 145-/xm quartz sample with reflections from each of the two surfaces. For 
this calculation, again Uc = 1.54661, rio = 1.53773, we ignore the frequency dependence of S, and the fast axis of the 
quartz plate is aligned with the horizontal axis in the laboratory frame so that ai = 0. The magnitude of the reflectance 
from each surface is assumed to be the same so that |rop = kiP- The SPDC spectrum is calculated explicitly via 
solutions to the phase-matching conditions using published Sellmeier equations for BBO. We are interested in the 
particular case of degenerate, coUinear type-II phase matching. The top two plots represent the expected rate of 
coincidence when the sample photon is mixed with a horizontally polarized (Rh) or a vertically polarized (Ry) 
reference photon. The bottom trace represents the re-normalized total coincidence rate [Rt — {Ry + Rh — Ao)/Ao] 
from which the relative reflectance and positions of each interface can be determined. 

In the Ry curve (middle trace), there is no dip at the flrst interface since the polarization mode reflected from 
this interface is solely horizontal. The polarization state is altered via propagation through the birefringent material 
and contains both vertically and horizontally polarized photons at reflection from the second interface. The peak 
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between the two interfaces in Ru (top trace) and Rt (bottom trace) is result of interference between each layer. This 
peak (which can alternatively become a dip depending on the phase accumulated between the layers) is susceptible 
to dispersion in the sample, unlike the dips that correspond to sample layers. Thus the dispersion properties of the 
material can be extracted from this feature. 

In summary, we ascertain that three experiments are required to completely determine the sample properties. 
We first select the reference arm polarization to be horizontal (H) and measure the quantum interferogram i?H by 
recording the coincidence rate of photons arriving at the two detectors as the path-length delay ct is scanned. The 
reference arm polarization is then rotated into the orthogonal vertical (V) polarization and a second measurement 
is made to measure the quantum interferogram Ry. The third measurement is made by selecting a value of cr that 
coincides with the position of a layer. The angles of the polarization elements in the reference arm are then adjusted 
to maximize the coincidence rate. 

The sample properties are found as follows: by forming a ratio of Ay and Ah at a value of cr that coincides 
with the position of a layer, we can determine the value of the birefringence contained in the parameter 6; using 
the angles from the polarization elements in the reference arm, a can be found from solving the equations for 
orthonormality. This technique is similar to nulling techniques in ellipsometry; and the total quantum interferogram 
Rt can be computed from the sum of Ru and Ry, then readjusted for the dc offset given by the constant term Aq, 
i.e. Rt = {Ry + Rh — Ao)/Ao. The Rt curve provides the path-length delay between the interfaces as well as the 
ratio of the relative reflectance from each layer. 

V. CONCLUSION 

We have set forth a new polarization-sensitive QOCT (PS-QOCT) scheme and provide a general Jones matrix 
theory for analyzing its operation. PS-QOCT provides a means for determining information about the optical path 
length between isotropic reflectors, the relative magnitude of the reflectance from each interface, the birefringence of 
the material between the interfaces, and the orientation of the optical axis a of the sample. Inasmuch as PS-QOCT 
is immune to sample dispersion, measurements are permitted at depths greater than those accessible via ordinary 
optical coherence tomography. 
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FIG. 1: Conceptual diagram of polarization-sensitive quantum-optical coherence tomography (PS-QOCT). The system is based 
on a Mach-Zehnder interferometer in which twin photons from SPDC, represented by the vector Jin, are separated into two 
arms at a polarizing beam splitter (PBS). The signal photon at angular frequency uj travels in the sample arm and experiences 
a path delay di as well as an arbitrary polarization rotation described by Ui. The reference arm contains the idler photon 
at angular frequency uj' which experiences a path delay d2 and a user-selected polarization rotation U2. Paths CD and @ 
impinge on a final beam splitter (BS) which mixes the spatial/polarization modes into paths (D and ®. Jout represents the final 
twin-photon Jones vector from which the fields at the detectors and the final coincidence rate are computed. 
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FIG. 2: Possible implementation of polarization-sensitive quantum-optical coherence tomography (PS-QOCT). A narrow-band 
pump laser at a wavelength of 400 nm pumps a 1.5-mm-thick /3-barium borate (BBO) nonlinear crystal (NLC) oriented for 
type-II, coUinear SPDC with a center wavelength of 800 nm. The pump beam is removed from the SDPC by use of a highly 
reflective mirror (HR 400) centered at the pump wavelength concatenated with a long pass filter (LP 695). The vertical and 
horizontal components in the SPDC beam are separated by a polarizing beam splitter (PBS) into the reference arm and sample 
arm of a Mach-Zehnder interferometer. The reference arm consists of a variable path-length delay comprised of a half-wave plate 
(HWP), a second polarizing beam splitter (PBS), a quarter- wave plate (Q), and a translational mirror. The final polarization 
of the reference beam (indicated as 0) can be oriented to either vertical or horizontal by a linear rotator prior to the final 
beam splitter (BS). The sample arm consists of a beam splitter and a quarter- wave plate (Q) so that circularly polarized light 
is normally incident on the sample. The back-refiected light from the sample, which in general has elliptical polarization, mixes 
with the delayed reference beam at the final beam splitter (BS). The outputs from the BS are directed to two single-photon 
counting detectors. The coincidence rate R{t) for photons arriving at the two detectors, as a function of the path-length delay 
CT, are recorded in a time window determined by a coincidence-counting detection circuit (indicated as (8)). 
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FIG. 3: Sample comprised of A'^ reflective layers. The probe beam is incident at the right. Each interface at position Zm is 
described by a reflectance matrix r,n(w). The optical properties of the sample layers between interfaces are described by the 
Jones matrix Sm- 
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FIG. 4: Simulation results for a single reflector buried beneath 120 fim of quartz with ric = 1.54661, Uo — 1.53773, and |rip — 1, 
using the scheme shown in Fig. |5] The optical axis of the quartz sample is aligned with the horizontal axis in the laboratory 
frame so that ai = 0. The top two curves represent the normalized coincidence rate when the sample photon is mixed with 
a vertically polarized (i?v, dash-dot curve) or a horizontally polarized {Rh, dashed curve) reference photon. The solid curve 
represents the total coincidence rate (-Rt) from which the reflectance of the layer can be recovered. 
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FIG. 5; Simulation results for a 145-/im quartz sample with reflections from each interface using the scheme shown in Fig.|5| For 
this calculation, Tie ~ 1.54661, rio — 1.53773, and |ro|^ — \rif ■ The optical axis of the quartz is aligned with the horizontal axis 
in the laboratory frame so that ai — 0. The top two plots represent the normalized coincidence rate when the sample photon 
is mixed with a horizontally polarized (-Rh) or a vertically polarized (-Rv) reference photon. The bottom trace represents the 
total coincidence rate (-Rt) from which the relative reflectance of each interface can be recovered. The separation of the dips 
is given by the optical path length of the quartz nL ~ 224 ^m. 



